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ABSTRACT 

We present Monte Carlo and Langevin micromagnetic calculations to investigate thermal 
switching of single-domain ferromagnetic particles. For the Monte Carlo study we place 
particular emphasis on the probability that the magnetization does not switch by time t, 
Pnot{t)- We find that Pnot(^) has different behaviors in different regimes of applied field, 
temperature, and system size, and we explain this in terms of different reversal mechanisms 
that dominate in the different regimes. In the micromagnetic study of an array of Ni pillars, 
we show that the reversal mode is an 'outside-in' mode starting at the perimeter of the 
array of pillars. 

INTRODUCTION 

All facets of the dynamics of nanoscale magnetic materials are currently active areas of 
research. The ability to construct single-domain nanoparticles via various methods and to 
measure the properties of individual nanoparticles and small arrays of nanoparticles [1 — 
4] provides clean experiments compared with previous studies of mixtures of magnetic par- 
ticles. One of the driving forces for applications is the rapid increase in density of magnetic 
recording devices, and the associated need to store each bit of information on a smaller 
number of grains [5]. On the theoretical and simulational side, a detailed understanding 
of nucleation and growth mechanisms that lead to the decay of a metastable state in finite 
systems has led to the identification of different decay modes in different parameter regimes 
[6 — 10]. In addition, new simulation algorithms are becoming available that may allow 
microscopic simulations at the inverse phonon frequency to extend to the technologically 
important time scales of years [11, 12]. 

In this brief paper we concentrate on using the thermal activation picture of nucleation 
and growth in simple metastable systems to better understand the reversal mechanisms 
for more realistic models of magnetism. We present Monte Carlo simulations for Pnot(^), 
the probability that the metastable magnetization has not yet reversed at time t. We 
show how well this simulation of a heterogeneous system fits our theoretical description for 
-Pnot(^)- We also present Langevin micromagnetic results for an array of Ni pillars. In the 
Langevin micromagnetic calculations the switching mechanism involves escape over a saddle 
point driven by random thermal fluctuations at constant field, rather than the deterministic 
disappearance of the metastable state as the field increases [13, 14]. Consequently, the 
Langevin micromagnetic calculations have the ability to measure Pnot(^) directly. 

MODELS AND METHODS 

We have preformed two types of computer calculations. The first consists of Monte Carlo 
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simulations [15] of the square- lattice Ising model with Hamiltonian 



where Si=±l and the local fields Hi are random numbers uniformly distributed between 
a maximum and a minimum value. This work uses periodic boundary conditions and a 
Glauber Monte Carlo update at randomly chosen sites. The initial state has all spins up, 
and at t—0 a negative field {Hi} is apphed. The unit of time is Monte Carlo Steps per 
Spin (MCSS). A rigorous derivation of the stochastic Glauber dynamic for Ising models 
from microscopic quantum Hamiltonians has been established under certain conditions in 
the thermodynamic limit [16], with the Monte Carlo time unit related to heat-bath phonon 
frequencies. This simulation has been performed to test our prediction for the forms of 
-Pnot(^) in a system with quenched bulk randomness. 

In order to simulate models with realistic spin degrees of freedom, we have programmed 
a Langevin micromagnetics code similar to that reported in [17]. We have used a phe- 
nomenological damping parameter a, and classical spins of constant length given by the 
bulk saturation magnetization Mg. Then at each lattice site i there is a scaled magnetiza- 
tion m — Ms/Mg. The standard Ginzburg-Landau-Lifshitz micromagnetic equation [13, 14] 
is 

The scaled effective field at each site, hi^^s, contains contributions from terms including 
the exchange interaction, the dipole-dipole interaction, the interaction due to crystalline 
anisotropy, the applied field, and a scaled noise term proportional to the the Langevin fields 
({t) [14, 17]. The Langevin noise term ^ and the integration time step At are related by 
(ocy^At. Even though the set of equations used in this Langevin micromagnetics simulation 
are approximations to the actual equations [18], the approximation should be reasonable far 
below the critical temperature. We have used a fourth-order Runge-Kutta algorithm as the 
integration scheme in order to keep the length of the rhi constant. 



RESULTS 

In the Ising simulation, we used 100x100 lattices at T=0.8Tc~1.815J. As in the case 
of homogeneous nucleation [7], homogeneous nucleation for single-domain Ising particles 
with demagnetizing fields [8], and for single-domain Ising particles with different boundary 
conditions [9], we have identified different decay regimes [6 — 10]. The functional form 
for quantities such as Pnot(^) are different in the different decay regimes. In particular, in 
the single-droplet (SD) regime of a single-domain magnetic particle, Pnot(^) = exp(— t/r). 
Here r is the average lifetime for the decay of the metastable magnetic state due to thermal 
fiuctuations. However, in the multi-droplet (MD) regime for single-domain particles the 
functional form for Pnot(^) is given by 

PnotW = ^erfc(^) , (3) 

where erfc is the complementary error function, and the width A depends on the system size. 
This form for P^ot (t) results from the assumption that many independent droplets nucleate 
and grow in different parts of the system, collectively leading to the magnetization reversal. 
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Figure 1: The probability Pnot(^) in the multidroplet (MD) regime for a 100x100 Ising 
model with periodic boundary conditions at T—O.STc. The data (open circles) are for 1000 
escapes from the metastable state. The random fields Hi are uniformly distributed, and 
centered on —0.34725 with a width of 0.34725. Two different one-parameter fits to the data 
are shown. The solid line is a fit to the complementary error function given by Eq. (3). 
The dashed curve is a one-parameter fit to a stretched exponential. In both of these fits the 
average lifetime was set at the measured value of t=54.4 MCSS. (So that the solid fine can 
be seen, only a small number of the 1000 data points are shown.) 

Figure 1 shows data obtained from 1000 reversals with the fields Hi uniformly distributed 
with width 0.34725 J centered about —0.34725 J. This distribution for Hi ensures that 
the system is in the multi-droplet regime [7]. This figure should be compared with Fig. 3c 
of Ref. [7] which is for the Metropolis dynamic with a uniform field. The hfetime r was 
measured to be t=54.4 MCSS, where the time unit is Monte Carlo Steps per Spin (MCSS). 
Using the Mathematica nonlinear fit function to fit Eq. (3) to the data, gives the value 
A=4.26 MCSS, which is shown as a solid curve in Fig. 1. It is also possible to try other 
standard expressions for the sigmoidal curve for Pnot(^) in the MD regime. Fig. 1 also shows 
a one-parameter fit to a stretched exponential, which gives Pnot(^) = exp (— 6t^°-^), shown 
as the dashed line. Here b is determined by requiring that the average of the distribution 
be r. Clearly among these one-parameter fits, the one to Eq. (3) fits the data much better 
than does a stretched exponential. 

Figure 2 shows an example of the type of thermal switching simulations [19] that can 
be performed using Langevin micromagnetic calculations. This figure represents a square 
array of magnetic Ni pillars. Similar arrays of Fe pillars have been built and measured 
experimentally [1, 2]. The simulation was started with all spins pointing up, and at t=0 
a field parallel to the spins was applied for 1 nsec to allow the system to come to thermal 
equilibrium. Then the field was reversed to point opposite to the average magnetization, 
leaving the spins in a metastable state. However, the magnitude of the applied field would 
not have been sufficient to switch the system at zero temperature, and the switching event 
depicted in Fig. 2 is enabled by the thermal fiuctuations included in Eq. (2). The random 
thermal field has its strongest effect at the edges of the array, where the demagnetizing field is 
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weakest. In particular, an 'outside-in' switching mode is seen, in which the decay towards the 
stable magnetizaton direction starts from the pillars at the edge and subsequently propagates 
to the pillars in the interior of the array. 

It is important to note that both the Monte Carlo and the Langevin micromagnetic 
simulations were conducted for applied fields sufficiently weak that a free-energy barrier 
against decay of the metastable state remained, and that this barrier had to be overcome 
by the thermal fiuctuations in order for the magnetization switching to occur. 

CONCLUSIONS 

Realistic simulations of models for magnetic materials have been carried out. The form 
for the probability of not switching, Pnot(^), is shown to fit the form of a complementary 
error function. The reversal mode for arrays of single-domain magnetic pillars has been 
identified as reversal of pillars on the boundary, followed by reversal of interior pillars. 

ACKNOWLEDGMENTS 

This research was supported in part by NSF Grant No. 9520325, FSU MARTECH, and 
FSU SCRI (DOE Contract No. DE-FC05-85ER25000). Supercomputer access was provided 
by the DOE at NERSC. 

REFERENCES 

1. A. D. Kent, S. von Molnar, S. Cider, and D. D. Awschalom, J. Appl. Phys. 76, 6656 
(1994). 

2. S. Wirth, M. Field, D. D. Awschalom, and S. von Molnar, Phys. Rev. B, Rapid 
Communications, in press (1998). 

3. W. Wernsdorfer et al., Phys. Rev. Lett. 77, 1873 (1996); Phys. Rev. B 55, 11552 
(1997); Phys. Rev. Lett. 78, 1791 (1997). 

4. M. Lederman, S. Schultz, M. Ozaki, Phys. Rev. Lett. 73, 1986 (1994). 

5. J. C. Mallinson, The Foundations of Magnetic Recording (Academic, New York, 1993), 
Second Edition. 

6. P. A. Rikvold, H. Tomita, S. Miyashita and S. W. Sides, Phys. Rev. E 49, 5080 (1994). 

7. H. L. Richards, et al., J. Magn. Magn. Mater. 150, 37 (1995). 

8. H. L. Richards, et al., Phys. Rev. B 54, 4113 (1996). 

9. H. L. Richards, et al, Phys. Rev. B 55, 11521 (1997). 

10. For a review see P. A. Rikvold, M. A. Novotny, M. Kolesik, and H. L. Richards, in 
Dynamical Properties of Unconventional Magnetic Systems , edited by A. T. Skjeltorp 
and D. Sherrington, NATO Science Series E: Apphed Sciences, Vol. 349 (Kluwer, 
Dordrecht, 1998). 

11. M. A. Novotny, Phys. Rev. Lett. 74, 1 (1995), Erratum 75, 1424 (1995); M. Kolesik, 
M. A. Novotny, and P. A. Rikvold, Phys. Rev. Lett. 80, 3384 (1998). 



4 



12. A. F. Voter, Phys. Rev. Lett. 78, 3908 (1997); J. Chem. Phys. 106, 4665 (1997). 

13. A. Aharoni, Introduction to the Theory of Ferromagnetism, (Clarendon Press, 
Oxford, 1996y 

14. W. F. Brown, Phys. Rev. 130, 1677 (1963). 

15. K. Binder, in Monte Carlo Methods in Statistical Physics, edited by K. Binder, 
(Springer, Berlin, 1979). 

16. P. A. Martin, J. Stat. Phys. 16, 149 (1977). 

17. E. D. Boerner and H. N. Bertram, IEEE Trans. Magn. 33, 3052 (1997). 

18. D. A. Garanin, Phys. Rev. B 55, 3050 (1997). 

19. G. Brown, M.A. Novotny, and P.A. Rikvold, in preparation. 

20. L. Neel, Ann. Phys., Paris 3, 137 (1948). 



5 



Figure 2: A scries of snapshots of a Langevin micromagnetic calculation for magnetization 
reversal in a square array of Ni pillars that are 200 nm tall, 200 nm apart, and each have 
a diameter of 40 nm. Each pillar is discretized using 5 lattice points. (For clarity of 
presentation, the vertical scale of this figure is enhanced compared with the horizontal 
scale.) The temperature is 300 K, the spins arc initially up, and the applied field is down 
with a magnitude of 1225 Oe. The integration time step is At=l pscc. The time sequence 
following field reversal, (reading line by line - left to right and top to bottom) is 4 nsec, 
8 nsec, 12 nsec, 16 nsec, 20 nsec, and 24 nsec. This reversal mode is different from coherent 
rotation [5, 20]. 
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